Disordered free fermions and the Cardy Ostlund fixed line at low temperature 
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Using functional RG, we reexamine the glass phase of the 2D random-field Sine Gordon model. 
It is described by a line of fixed points (FP) with a super-roughening amplitude (w(0) — u(r)) 2 ~ 
A(T) In 2 r as temperature T is varied. A speculation is that this line is identical to the one found in 
disordered free-fermion models via exact results from "nearly conformal" field theory. This however 
predicts A(T = 0) = 0, contradicting numerics. We point out that this result may be related to 
failure of dimensional reduction, and that a functional RG method incorporating higher harmonics 
and non-analytic operators predicts a non-zero A(T — 0) which compares reasonably with numerics. 

PACS numbers: 75.10.Nr, 72.15.Rn, 64.60.Ak 



I. INTRODUCTION 



A. Overview 



Few analytical methods can handle systems with 
quenched disorder beyond mean field and phenomeno- 
logical arguments, and each seems adapted only to a 
specific class of models. One example is the functional 
RG (FRG)i used to study elastic objects in a random 
potential 2 ^ or continuous (e.g. XY) spin models in a 
random field^. When pure, these systems possess internal 
order and one studies the case where the (weak) quenched 
disorder produces only elastic deformations (no topolog- 
ical defects) parameterized by a displacement field u(x) 
away from perfect order u = 0. They exhibit glass 
physics and pinning, i.e. thermal (and quantum) fluctu- 
ations are small and energy minimization yields multiple 
metastable states. The FRG handles this by summing an 
infinite set of operators, all relevant for internal dimen- 
sion d < d uc , usually d uc — 4. The coupling constant is 
a function, R(u), which represents the correlator of the 
coarse-grained disorder. FRG fixed points R*(u) were 
obtained in a d = d uc — e expansion for various disor- 
der classes^ and recently it was shown how to directly 
measure them in numerics 5 . 

Disordered systems in d = 2 can sometimes be stud- 
ied using methods inspired by conformal field theory 
(CFT). Exact results were obtained for free fermions 
with disorder of various symmetries using bozonisation 
and supersymmetry®' . The physics there is localiza- 
tion, and in recently much studied symmetry classes, 
de-localization near a point where the density of state 
(DOS) is singular (e.g. bipartite lattices leading to Dirac 
fermions). There, sample (and local DOS) fluctuations 
can become broad leading to freezing transitions in the 
DOS dynamical exponent z, p(E) ~ E~ 1+2 / z , as found 
for Dirac fermions with random vector potential^. Be- 
ing equivalent to a (dilute) limit of a random gauge 



XY spin model it exhibits an infinite set of relevant 
operator a 9 ! 10 which can be handled using another func- 
tional RG method 11 based on the KPP-Fisher equation 
(hence noted here KPP-FRG). Freezing transitions are 
expected to impact a wider class of models, e.g. the 
time-reversal, particle-hole symmetric class, where the 
p(E) ~ E~ 1 e~ c ( lnE ' > ' i Gade-Wegner singularity (p = 
1/2) was argued to become fi = 2/3^2. This was 
first found using KPP-FRG on the dual spin version 
(the Cardy-Ostlund model with dilute vortices)^ and, 
recently, within the fermion model 7 ^ and the sigma 
models. These emerging 2d field theories still need test- 
ing on whether they capture all (nonperturbative) freez- 
ing (and glassy) physics, for which the spin models can 
provide insight. 



B. Model 

In this paper we focus on the simpler 2D random field 
XY (Cardy-Ostlund) model without vortices^ . Being 
the bosonized form of a fermion model, exact results were 
claimed for this model using CFT inspired methods^. 
Since it is also a (periodic) pinning problem, we can also 
use FRG and test the strong disorder regime, expected to 
be correctly described already to one loop. The model, 
also called random phase Sine Gordon model (refered to 
here simply as CO model) also describes a periodic object 
(such as an array of flux lines) at equilibrium on a 2d dis- 
ordered substrat o 2 ' 17 at temperature T, and of partition 
function Z = J Vue~ Sco ^ with ScoM = H co [u]/T 
and: 



H C o[u] = j d 2 r[|(V r u) 2 - h ■ V ' r u - Re(£e 



')] (1) 



where u(r) G] — 00, +00 [ is an XY phase (excluded vor- 
tices), k the elastic coefficient and £(r) a (complex) 
quenched random field with random phase £(?")£* (r') = 



2 



2g S 2 (r — r'). The additional random field h(r) is gener- where f = J dx J Q dr, ip = iJjl + ipR, P = ip^ip i 
ated under coarse-graining with h l (r)hi (r') = o~8ij5 2 (r — 
r'). Disorder being short ranged, the model ([lj, to- 
gether with some of its discrete versions, exhibits the 
STS symmetr y 19 ' 20 , i.e. statistical invariance of the non- 
linear part under u{r) — ► u(r) + v(r). The model ([1]) 



the fermion density and g 2 ,4 are (spinless) fermion in- 
teraction parameters with p s = iplip s , s — L,R. Using 
bosonization: 



exhibits a super-rough phase for T < T g = inn with 
anomalous growth of the 2-point correlatio n 16 ' 18 ' 19 : 



B(r) = ((u(r) - u(0)) 2 ) = A(T)(\nr) 2 + O(lnr) (2) 

where (...) denotes thermal averages. This is in con- 
trast with the usual logarithmic roughness for T > T g , 
disorder being irrelevant, i.e. A(T > T g ) = 0. The 
STS symmetry implies that T/T g is uncorrected by dis- 
order and does not flow under RG. This "marginal glass" 
below T g is thus described at large scale by a line a 
fixed points (FP), indexed by T/T g (which can be iden- 
tified from the large r thermal connected correlation 
((u(r) - w(0)) 2 ) c - 4T/T g \nr). The line of FP can 
be reached near T g by perturbative one loop RG which 
predicts^: 



A(T) 



2r 2 



for small 



r={T g - T)/T g « 1 
(3) 

in agreement with precise numerics^ 2 -. At low tem- 
perature, only numerical results are available at 
present: it was concluded from exact ground state 
determinations 2 ^! for a closely related solid on solid 
model that super-roughening holds at T — 0, i.e. A(T = 
0) « .5 > (see below). 



II. RELATIONS TO DISORDERED FERMIONS 

A. General framework and bosonization 

A recent study 6 - of a free fermion model opened hope 
that A(T) could be obtained exactly. It has been known 
that weak disorder on the CO model at any T can be 
equivalently studied by adding quenched (dynamical) dis- 
order to a (euclidean) d = 1 + 1 interacting fermion 
model 2 ^ . Starting with the CO model without disor- 
der, i.e. h,£ = in |T]), one defines r = (x,t) where 
t G [0,/3ft] is imaginary time, and u(r) = 2<p{r). Then 
the bosonic action Sb [4>] — 6>co [ u — 20] , where one sub- 
stitutes y V 2 — > 2^f(«^t + v ^l) m ©> is tne euclidean 
action (with Z = J T>4>e~ Sb ^) of the bozonized form of 

the fermion model (with Z = j Vip^Vipe^f^^^) of eu- 
clidean action Sf = So + Smt with: 

So = f bP R d T i>R + + ~ f Ho (4) 



H = hv F : ip[(id x - k F )ipL - ip R {id x + k F )i/j R : 

J X 

Sint = ^ J y (arAr + PlPl) + 92PrPl (5) 



V27ra 

with s = 1/ — 1 for R/L movers, the density reads: 

p = po - d J± + J_Re(e 2 ^-~ 2 ^) 
7r 2na 



(6) 



(7) 



where po is the average density and a a UV cutoff. The 
(Tomonaga Luttinger) model So + Sint is equivalent, us- 
ing 27rp s = —d x 4> + sd x 9, to the quadratic action Sb[4>] 
and to the CO model without disorder, with: 



K = T/T g 



(8) 



hence any temperature can be reached by varying the 
Luttinger parameter K, i.e. the interactions, with: 

K 2 = (1 + 2/4 -2/ 2 )/(l + 2/4 + 2/2) (9) 



and the sound velocity 



v F 



1 



2/4) 



2/1 



(10) 



(11) 



which can be set to one by a rescaling of space-time. 
Next, the (formal) perturbation of the fermion model 
Sf — > Sf + Sdis with 



Sdi 



Wp 



W(r) = /*(r) + f(r) 



2'ikr. 



+ f(r)e 



— 2ikpx 



(12) 
(13) 



by a small (dynamic) random potential W(r) (neglect- 
ing for now higher harmonics) is equivalent via Sb[4>] = 
Sco[ u = 20], to the disorder pertubation in the CO 
model, i.e h,£ in ([1]). The correspondance reads h x = 
—p/(2nh), h T = 0, y£( r ) = £(r)/(wafi). An isotropic 
h distribution is obtained by coupling a random field to 
d T 4> but it makes no difference at large scale. For T < T g 
small disorder go is relevant and grows under RG to finite 
value on the fixed line 25 . At a naive level (see below) the 
relation © still holds since being unrenormalized (for 
exact STS). 



B. Predictions from "nearly conformal" field 
theory 

The recent study 6 - adresses a free fermion model in 
presence of (isotropic) disorder Sff = So + S' dis . We 
follow the notations of Ref 7 which uses disorder param- 
eters gA — tto-/T 2 and g&i cx go, shifted fermion fields 
■0s = e l ( e ~ s ^ and the same bosonic fields 4>l/r — 4* — s6. 
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One has chosen H = 1 = Vp , h x = A t /tt, h T = —A x /tt. 
This model was also found to be described by a line of 
FP with the exact result (denoting £ = l/(l + <?jvf/7r) and 
4> = u/2 the same bosonic phase as above): 



dign = , dig A 
(4> a (z,z)ct> b (0)) = 



(C9m) 2 



2tt 2 



9 2 mC 4 
{2-Kf 



In 2 zz 



(14) 

C,{5ab + — 
7T 



for real copies a, 6 of the model (no replica). The vanish- 
ing beta function for gu strongly suggests that the free 
fcrmion model is exactly on the CO fixed line. From the 
above it means that the value of K is shifted from K = 1 
(no disorder) down to K — 1/(1 + g^ /it) (in presence 
of disorder), to all orders in gu according to Refi^i To 
lowest order in <?m one easily finds such a shift on model 
5*0 + Sdis, e.g. considering the replicated action 



So[i> a ] ~ / (9mO m 

r.t „7,_ „7,t 



9aO / 



lib 



(15) 
(16) 



While bosonization of the terms a ^ b yields the 
cos(2(0 a — <p b )) disorder terms of the replicated CO ac- 
tion (see[23lbelow). the a = b term yields an interaction^ 
as in Eq. (TJJ with 172 = 2gM, compatible to lowest or- 
der with K = 1 — (72/71" • To understand higher orders one 
notes that the interaction used in&£ is fully isotropic, and 
does not correct v (while gi corrects v to order g 2 ), hence 
involves also a 2k f coupling to current. Note that density 
only couplings preserve vK to all orders. In bosonized 
form it reads 



9M 

2tt 2 



J aaJ aa 



9m 
2tt 2 



Pt0) 2 + {d x cj>f 



(17) 



where J ab = i>Lb^t a - 

The bosonised theory of Ref£ hence exhibits super- 
oughening, and is "nearly conformal" as the h random 
field (the gA sector) is a (Larkin type) random force triv- 
ially eliminated by a shift using STS. The S ab term in (TT4^) 
gives thermal connected correlations hence one identifies 
C = T/T g . It implies the prediction (with r = 1 - T/T g ): 



A ff (T)=2r 2 (l-r) 2 



(18) 



consistent with one loop at small r, exhibiting a max- 
imum at t = 1/2 where Aff = 1/8 (see Fig. [T]) and 
vanishing at zero temperature. Clearly Aff(T = 0) = 
is incompatible with numerical results and equivalence of 
the models. The reentrant shape of the Aff(T) is sur- 
prising as one expects stronger effect of disorder at low T. 
This is reminiscent of reentrant phase boundaries found 
in the random gauge XY model, due to neglecting op- 
erators (moments of vortex fugacity) which become rele- 
vant at low T— Here higher harmonics of the exp(i0) 
field are known to be generated as an interplay of inter- 
actions and disorder. This can be seen in the Haldane 



o 
11 




FIG. 1: Right y-axis: A ff (T) given by formula (|T8j) ob- 
tained, as explained in the text, from the theory of Ref.— and 
yielding Afj(T = 0) =0. Left j/-axis: the exact ground state 
determination of disordered solid-on-solid models yield A(T = 
0) = a 2 = 0.57 from Refi^ and A(T = 0) = 2(2tt) 2 B = 0.51 
from Ref.—, shown together with the value obtained by 
the present FRG method. 



representation^: 

P(r) = (Po - -<9z<Mr)) V. 

7T ' ■ 



,2ipp x -2iptj>(r) 



(19) 



which yields higher harmonics cos(2p((f> a — <fi b )) in (|23p . 
Note that the isotropic disorder J pW + jV where j 
is the exact current, d T p + d x j = 0, produces a seem- 
ingly irrelevant coupling d x 4>(r)Re(e^ 2l ^ r ' 1 £(r)) . Upon 
replicating its aa contribution generates a shift in K, 
and none in v. Whether the higher harmonics are gen- 
erated in a free fermion model (e.g. as in*) properly 
taking microscopic cutoff into account deserves further 
investigations^. Weak disorder counting shows them to 
be irrelevant for 1/2 < T/T g < 1 but these dimensions 
could change along the fixed line well before that. We 
now study the effect of higher harmonics using the FRG. 



III. FUNCTIONAL RG 



A. Effective action and general structure of the 
RG flow 



We start by noting that Aff(T = 0) = is reminis- 
cent of a dimensional reduction result, i.e. an artefact 
of considering an analytic disorder R(u). This suggests a 
reexamination of the model at low T using FRG. We now 
indeed show that the cusp in R"(u) is necessary to get 
A(T = 0) > 0. Using replicas u a (r) = 11%, a = 1, .,n and 
averaging over £, h we obtain the bare replicated action 



4 



(with Z n = J Vue~ So ^): 

H ro Ju\ K 



So[u] 



T 2T 



E / Ku-MT 1 (2°) 



2^E/ i^O ^ lv< b Vu?G (u?)] (21) 
= < - u», / p = Jd d r, f = fd d q/(2n) d . The 



2 c(q 2 /2Af,q 2 /2Al) where A 



with = w, a 
bare propagator is 
and A; = Aoe - ' are UV and IR cutoffs and c(z, s) is a 
smooth cut-off function. We denote c(z, oo) = 1 — c(z) 
with c(0) = 1 and c(oo) = 0. At the bare level, Rq(u) = 
go cos(m) and Go(u) = <7 in (f2"5)) . The effective action 
Ti[u] (i.e. the Legendre tranform of the free energy in 
presence of sources) admits a double expansion in number 
of replica (cumulants) and gradients, which starts as 



i l J 2T ^ 



x u - q {yi) 



(22) 



as in (|20I23[) with (periodic) renormalizcd function i?; (u) 
and Gi(u). Note that the STS symmetry constrains the 
form of the one replica term to be identical to the bare 
one. In d = 2 a p-periodic Ri becomes relevant below 
T = T 3 (p) = T g 2Tr/p and at T = is relevant below 
d — 4. Higher cumulants are irrelevant compared to Ri 
near d = A. In d = 2 a 1-loop approximation keeping 
only Ri yields the exact and reasonable properties 
at low T 31 i 32 . A specificity of d = 2, not systematically 
studied with the FRG (see however—) is that the term 
Gi must be retained. Indeed at T = and d — 2 the 
whole function Gi(u) is dimensionless at the bare level, 
and acquires anomalous dimensions due to temperature 
and disorder. The 2-point correlation function reads: 



T 



(u q U^ q ) = ^-Gf 



^!(-^'(0) + G / (0)g 2 ) (23) 



Near T g , it was shown that Gi (0) ~ I which yields, setting 
I = log(l/g), the In 2 |r| behavior in B(r) ([2]). We now 
investigate how a similar behaviour, G;(0) oc / can arise 
also at low T where, as for Ri(u), all harmonics in Gi(u) 
must be kept. We first set k = 1 in Eq. ([2H1) . Defining: 



T 



S d A?~ 2 T 



Rl(u) = S d Af- A Ri(u) 
G(u) = S d Af- 2 Gi(u) 



(24) 
(25) 
(26) 



with S2 — (2?r) 1 , the coupled RG equations for Ri and 
Gi have the structure 



d l T = {e-2)T , d l R{u)=f3 k {u) 
dtGiu) = (e-2)G(u)+fls(u) 



(27) 



with e = 4- d, fag = f3 R6 [R,R , R'" , G, G'...] where 
R"(tt) = If" (it) - lf"(0) and G(u) = G(u) - G(0). Under 
the following assumptions : i) G(m) is continuous in u = 
0, ii) Pq(u) has a good limit when u — > + : 



^(0+) = A d 



(28) 



and iii) G(u) has a good fixed point for all d one obtains 
that G(0) = Ad/ {2 — e), where admits an e-expansion 
around d = 4. Thus G(0) is diverging for e — > 2, yielding 
the announced behavior in d — 2, G(0) = A2I. Note that 
this implies: 

j^„m„+MM (29) 

y °2 y 

hence, going back to real space A(T = 0) = A 2 . 

B. Lowest order calculation 

We first compute, at T = 0, /3a, #5 to 0(R 2 ). More 
details are given in Appendix [X] From one loop 1PI 
graphs one finds = Sq[u] + ST[u]: 



ST[u] 



AT 2 



IE / (5r V ) 2 K« b )K«) (30) 

ab Jx V 



excluding terms proportional to T (tadpoles and some 
3- replica terms). Differentiating the local component of 
5T[u] in (|30| w.r.t. the IR cutoff A/ at fixed i?o, and 
inverting Rq = R t + O(Rf), one obtains the standard 
T = 1-loop FRG equation for Ri in d = 4 — e: 

= eiZ,(u) + (1/2)(R») 2 + 0(i?G, G 2 ) (31) 

with Ri(u) = 2I qRi(u) where 

I nm = S^Ai d+i+n - m ( y^Gf^g? (32) 
Jy 

for n,m even (see also Appendix |B|) . Eq. (j3"Tj) admits a 
non-analytic, p-periodic fixed point solution 2 - given by 

R*( P u) = (ep 4 /(72))(l/36 - u 2 (l - uf) (33) 

with u € [0, 1[, yielding the amplitude of the cusp: 

R*"'{p + ) = -R*"'(0-)=pe/G (34) 

Next, the bilocal component of ST[u] in (I30j) . when ex- 
panded in gradients generates a term VM° b Vw° b G;(M° b ). 
Proceeding as above one obtains 3 ^ (see Appendix up 
to C(i?G,G 2 ): 

diGiiu) = (e - 2)G(u) + (l2o/d)R"'(u) 2 (35) 

Eq. (05]) yields, in d = 2, G(0) - (I 2 o/2)(-R*'"(0)) 2 Z 
which has a non ambiguous value (f3~4l . It also leads to: 



G(«) = (/ 20 /d)((i?'"( U )) 2 - (i?'"(0)) 2 )/(2 - e) (36) 
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thus satisfying the requirement iii) above for e < 2. As 
shown below, the divergence for e = 2 is cured by consid- 
ering higher order terms in Eq. (|35| . From (f23l) . (f35l> a 
first estimate A {1) of A(T = 0) in © up to order 0(R 2 ) 



aw = ^(ir"'(o+)) a ' r '~ h " 



288 J 00 



(37) 



As announced it is non zero at T — only because the 
fixed point (|34p is non- analytic. The p 2 dependence in 
Eq. Q37p yields that is independent of k. Indeed 
one has B K , p (r) = n~ 1 B l p /^{r) if B KtP (r) denotes the 
correlation in ([2} for model ([1]) with p-periodic potential 
(and h = 0). While A<2=2 is expected to be universal, 
our truncation (|37|) depends on the shape of the IR reg- 
ulator c(z), with local sensitivity minimum around the 
exponential, c(x) — e~ x which yields A^ = 0.881. 



C. Higher order calculation 

We now push the T = analysis further by computing 
PfcPa up to order 0(i?G,G 2 ), which are of order 0(1) 
in d = 2. From ST[u] to one loop, one finds: 



0T=° = eR + I 00 R" 2 + 2/ 02 R"G + Io 4 G 2 ... (38) 



R 

G =Q = {ho/d)R" r * + 7l i?'"G' + 2/ 00 R"G" (39) 
+2/ 02 GG" +72G' 2 + ... 

where (see Appendix [Bl for more details): 

ll =4(I oa +K Q /d + I 22 /2d) (40) 
72 = 6/02 + I 2i /d + K 2 /d (41) 

S^A- d+i - n I y l d l Gf\7 J \7 n gf = SijKn/d (42) 

J V 



Remarkably, (j3"5l |3"9")) admit a simple solution: 

R"(pit) = p 2 ad,u(l — u), G(pu) = p 2 aadu{\ — u) (43) 
with: 

(44) 



a-d 



12(I oa + 2aI 02 +a 2 I 04 ) 



and e = 4 — d. Eqs. (|3"9"|) . (|4"3")l satisfy all above require- 
ments i), ii) iii). It yields a value for A4, up to order 
0(RG, G 2 ) at one loop, for a p-periodic potential: 



A d 



P e -W^ + «7i + Q 72 
144 (loo + 2aio2 + a 2 / 4) : 



(45) 



yielding back (|37[) if a is (formally) set to zero. Here a is 
the solution of a cubic equation for e < 2 (with a = 0(e) 
near d = 4), and of a quadratic equation in d = 2 



a 2 (io2 + 72 ) + a(J o + 71) + = 



(46) 



Among the two solutions of Eq. (|46]) . only one is physical, 
and yields G(0) cx A 2 l and an estimate A (2) of A(T = 0) 
up to order 0(i?G, G 2 ) at one loop. Using relations valid 
in d = 2 it can also be written as: 



2 2 

pea 



02 



no 



144 {I m + 2aI 02 + a 2 I M f 



(47) 



A^ is independent of k and rather stable as the IR cut- 
off c(z) is varied. A Gaussian c'(x) cx e - ^ -1 / 2 ) ^ 2a > 
reaches a local maximum with a = 1/4 and A^ = 0.29. 
Thus the numerical results, A(T = 0) = a 2 = 0.57 of 
Refi2£ (notice that 02 refers here to the notation used 
in Ref. ? ~ and is different from 0^=2 as in our Eq. (|4"4")l ) 
and j4(T = 0) = 2(2it) 2 B = 0.51 of Ref.24 seem to lie in 
between our estimates A^,A^' (see Fig. Q]). 



D. Extension to low temperature 

To be convincing the scenario should be stable to a 
small temperature T. To one loop (|2"7| become: 

diR(u) = fi?"(u) - J 2 TG(u) + /3| =0 (u) (48) 
9;G(u) = (e - 2)G(«) + TG"(u) + /3g =0 (u) (49) 

where J 2 = (S^)- 1 Af~ 4 <9 ; V 2 £f =0 . At low T, and dimen- 
sion d = 2, we look for a solution of Eq. (|4"9")) where for 
it 3> T the functions i? and G are given by the T = 
solution (|4"3"| . while for small u ~ T there is a thermal 
boundary layer (TBL) solution of the form 

R"(u) = Tr"{u/T) , G(u) = A 2 Z + c + Ts(u/T)(50) 

where, with no loss of generality, g(0) = r"(0) = 0. The 
functions r"(z) and g(z) satisfy the coupled equations 
(for p = 1) 



= -a 2 —+r"- J 2 g + I 00 r" 2 + 2I 02 r" g + I 04 g 2 
6 



A 2 = g" + ^§V" 2 + 7l r'V + 2/ 00 r".g" 
a 

+ 2/ 2.g.g" + 72ff' 2 



(51) 



where we have used R"(0) = —ad/Qp 2 and the notations 
r = r(z),g = g(z). From Eq. (|5T]) . one easily checks the 
large \z\ behaviors 

r"{z)^a d \z\+r + O(\z\- 1 )) (52) 
g(z) ~ aa d \z\ + g + OQzl- 1 )) (53) 



which match correctly the cusp (|43|) outside the TBL. 
The constants are ro = ( J2/02 + /o4)/(2/ 2 2 — 2/o4^oo) and 
30 = -(l + 2r / o)/(2/ 2). 

We have solved numerically these RG equations (fSTj) 
and checked that it admits a well defined solution : this is 
depicted in Fig. [3]where we plot r"(z)/a 2 \z\ (respectively 
g(z)/aa 2 \z\) as a function of z in Fig [3] a) (respectively 
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z 



IV. CONCLUSION 

To conclude we applied the Functional RG to the ran- 
dom phase Sine-Gordon model, when gradient expansion 
is important. We showed that the super- roughening am- 
plitude A(T) does not vanish at T ~ thanks to the non- 
analyticity of the fixed point function i?(u), which avoids 
dimensional reduction. The one loop estimates give an 
order of magnitude consistent with numerics. Our study 
suggests that either the mapping between free fermions 
model and the CO model fails below some temperature 
or that new operators become relevant, e.g. higher har- 
monics, as is the case for interacting fermions models. 
Such a scenario was demonstrated in random field spin 
models^. Whether the path integral transformations of 
Ref& hold at low T needs to be clarified. A procedure 
exists 5 to measure the FRG functions described here, in 
numerics and check the scenario. We hope that it allows 
progress in issues common to fcrmion models, localiza- 
tion and periodic pinned systems. 
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FIG. 2: Numerical solution of the coupled RG equations (|5ip 
in dimension d = 2, for p = 1. We have used an exponential 
cutoff function c(x) = e~ ax . a) : Plot of r" {z) I ai\z\ as a 
function of z > 0. b) Plot of g(z)/a.a,2\z\ as a function of z. 
Notice that the convergence towards the asymptotic behavior 
(|53l) is quite slow. We have checked that both curves converge 
to unity at very large \z\. 

Figf2]b)). Hence this shows that the scenario is consis- 
tent: in fact these one loop equations, if naively continued 
already predict a In 2 behavior Q at all T and the cor- 
rect T g , with, near T g , R(u) ~ cosu since only the lowest 
harmonic is relevant. A more complete interpolation of 
A(T) deserves further investigations. 



APPENDIX A: GRADIENT EXPANSION 

The calculation of of the effective action up to order 
i?p yields: 

8T[u] = —[J f (M< b + 5< b )) (Al) 

ab Jx 

+ 4^2 E/ {R*{uf + 8uf)R {uf + 8uf))] 

abcd Jx V 

where (...) means average over the Su with the free prop- 
agator, keeping only the 1PI graphs. This gives: 



<ft>] = / R o« b )^l + I (er v ) a E 5 R 0« 6 )Ro«) + E <«X'K1 ( A 2) 

L± ab Jx Jx V ab L abc 

I 

where Rq(u) = Rq(u) — i?o(0). The last term is a third Consider a general bilocal functional: 
cumulant proportional to temperature, and here we ig- 
nore it. We now expand the second term in gradients. 

/ F(u x ,u v ,y ~ x) (A3) 

J xy 
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We can either use z = (x + y)/2, t = y — x and expand 
everything to second order in t : 



izt 



F{uz + \{t ■ V)w, + i(t • \7) 2 u z ,u z - l(t ■ VK + i(t • V) 2 u z , t) 



\{t ■ V)u*(t • V)u z (5i - a 2 ) 2 F(ii z ,u z , t) + i(i ■ V)(t ■ V)««(9i + d 2 )F{u Zl u z , t) 



(A4) 
(A5) 



Using also that J u, t) = (we assume that the local 
part has already been substracted) and performing inte- 
gration by part to get rid of the VV«F(it, u, t) term, we 
obtain a term containing (d\ — d2) 2 — {di+d2) 2 = — 4<9i<9 2 
and thus finally one obtains the gradient expansion as: 



/ 

Jxy 



F(u x ,u y ,y - x) 



= - [ \(t-V)u z (f^)uzdid 2 F{u z ,u z ,t) (A6) 
Alternatively one can also perform the expansion as 

/ F(u x ,u y ,y - x) (A7) 

Jxt 



= I F{u x ,u x + (t-V)u x + ht-Vfu x ,t) 

= I \[{t-V)u x {t-w)u x d 2 2 

Jxt 1 

+{t-V){t-V)uxd 2 ]F(ux,Ux,t) 



(A8) 



and integration by parts then yields the same result as 
above. Thus one finds that the part of interest of the 
effective action is: 

5I >] = ^E / <K b )KK b ) 

ab Jx 

x / (G 2 t 5(t) I g 2 )(t ■ VK 6 (t • V)< b (A9) 

Jt J s 

yielding the result given in the text. 



APPENDIX B: CALCULATION OF INTEGRALS 

In this appendix, we give the detailed calculation of the 
coefficients I n , m ,K n entering the RG equations given in 
the text. One introduces the notation Qf = Af _2 7, Aiy . 
Focusing on the large I limit, one defines = g v and 

dg v as: 

g v = f e <™_L(i _ c(g 2/ 2 )) , dg* = - [ c'(q 2 /2)e^ 

Jq Q Jy 

(Bl) 

Let us start with I ao : 

J- f^VnV = - 1 f C V/2), 

s d 



In a similar way, one computes the coefficients (in d = 2) 
I02 = -\ (B3) 

/>0O 

I Q4 = -2 / duuc'(u)(l - c(u)) (B4) 
Jo 

Next, we compute 7 2 o 

I*> = i- f y 2d 9 v 9 v (B5) 

&d Jy 

" |/r'V/2)(l-o( 9 V 2 )) 
+^ /c"'(9 2 /2)(l-* 2 /2)) 

/>oo 

= / du(l-c(«))(c"(«)/« + c"'(u)) in d = 2 
Jo 

Similarly, one has (in d = 2) 

/•oo 

/ 22 = -2 / du(l-c(u))(c"(u) + uc"'(u)) (B6) 
Jo 

/•oo 

7 24 = 4 / du(l-c(u))(ut/'(u)+uV'(u)) (B7) 
Jo 

Finally one has 

K = ^ [ d g y(vi g y)yi 

&d Jy 

= -j-J^'( q 2 /2)(l-c( q 2 /2)) (B8) 

= / due" {u){l - c{u))) in d = 2 (B9) 
Jo 

and similarly, in d = 2 

/•OO 

K 2 = 8 duuc" {u){\ - c{u)) (BIO) 
Jo 



- du— i-i(l- c («)) in d = 2 (B2) 
2 Jo u 
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